The Surface Tension of Quark Matter in a Geometrical Approach 



Marcus B. Pinto, 1 Volker Koch, 2 and J0rgen Randrup 2 

Departamento de Fisica, Universidade Federal de Santa Catarina, 88040-900 Florianopolis, Santa Catarina, Brazil 
2 Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 

(Dated: July 24, 2012) 

The surface tension of quark matter plays a crucial role for the possibility of quark matter nucle- 
ation during the formation of compact stellar objects, because it determines the nucleation rate and 
the associated critical size. However, this quantity is not well known and the theoretical estimates 
fall within a wide range, 70 ~ 5 — 300 MeV/fm 2 . We show here that once the equation of state is 
available one may use a geometrical approach to obtain a numerical value for the surface tension 
that is consistent with the model approximations adopted. We illustrate this method within the 
two-flavor linear a model and the Nambu-Jona-Lasinio model with two and three flavors. Treating 
these models in the mean-field approximation, we find 70 ~ 7 — 30 MeV/fm 2 . Such a relatively small 
surface tension would favor the formation of quark stars and may thus have significant astrophysical 
implications. We also investigate how the surface tension decreases towards zero as the temperature 
is raised from zero to its critical value. 

PACS numbers: ll.10.Wx, 12.39.Ki, 21.65.Qr, 26.60.Kp 



I. INTRODUCTION 

Lattice-gauge calculations yield a non- vanishing value of the quark condensate (ipip) in the QCD vacuum 
indicating that chiral symmetry is broken. This general feature of the vacuum remains present even for massless 
quarks because the symmetry is then broken spontaneously. On the other hand, chiral symmetry is expected to 
become restored at sufficiently high values of the net-baryon density p or/and the temperature T. The character 
of this phase change is not yet well understood but it has significant implications in areas such as cosmology and 
astrophysics and it is a focal point for current experimental and theoretical research in nuclear physics. 

Nuclear collision experiments carried out with the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven 
National Laboratory and with the Large Hadron Collider at CERN explore systems having relatively small net- 
baryon densities p and the associated chemical potentials p are negligible. Lattice calculations can readily be carried 
out at vanishing p and they indicate that a cross-over transformation from the chirally broken phase to the restored 
phase occurs as the temperature is increased from below to above the cross-over temperature T x rj 160 MeV [IHl]- 

The other extreme region of the QCD phase diagram, namely low temperatures and high chemical potentials, cannot 
be addressed by current lattice-QCD methods, due to the fermion sign problem, and studies of this phase region must 
therefore rely on less fundamental models. Most investigations suggest that there is a first-order phase transition 
which, for T « 0, sets in at baryon densities several times that of the nuclear saturation density, po ~ 0.153/fm 3 . The 
properties of strongly interacting matter in this phase region are important for our understanding of compact stars. 

If indeed such a first-order phase transition exists at T = 0, then, as the temperature is raised, one would expect it to 
remain present but gradually weaken and eventually terminate at a critical point (p c , T c ). The existence and location 
of such a critical point is a subject of intense theoretical investigation with a variety of models, including in particular 
effective-field models, such as the linear a model (LSM), and effective quark models, such as the Nambu-Jona-Lasinio 
(NJL), at different levels of sophistication considering up to three quark flavors and possibly including the Polyakov 
loop to account for confinement 0, [E| • Experimentally the corresponding region of density and excitation may be 
produced in current nuclear collisions at the low-energy end of RHIC and in the future with FAIR at GSI and NICA 
at JINR which are being constructed with such investigations in mind. 

In the present work, we concentrate on the high-/i and low-T part of the phase diagram with the aim of exploring 
the expected chiral phase transition which has significant implications for the possible existence of quark stars [fj|, l7| . 
It should be noted that chiral symmetry may be restored already during the early post-bounce accretion stage of a 
core-collapse supernova event and the associated neutrino burst might then provide a spectacular signature for the 
presence of quark matter inside compact stars 0. However, as pointed out in Refs. d, [l(|, the possibilities depend 
on the dynamics of the phase conversion and especially on the time scales involved. 

When the phase diagram of bulk matter exhibits a first-order phase transition, the two phases may coexist in mutual 
thermodynamic equilibrium and, consequently, when brought into physical contact a mechanically stable interface 
will develop between them. The associated interface tension 77- (which we shall often refer to simply as the surface 
tension of quark matter) depends on the temperature T; it has its largest magnitude at T — and approaches zero 
as T is increased to T c . This quantity plays a key role in the phase conversion process and it is related to various 
characteristic quantities such as the nucleation rate, the critical bubble radius, and the favored scale of the blobs 
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generated by the spinodal instabilities [ll|, ll2j. (As we shall see, the surface tension is essentially proportional to the 
effective interaction range, which determines the width of the surface region, and the spatial size of the most rapidly 
amplified density irregularity is also proportional to this quantity.) 

Unfortunately, despite its central importance, the surface tension of quark matter is rather poorly known. Estimates 
in the literature fall within a wide range, typically 70 « 10-50 MeV/fm 2 Q1GJ] and values of 70 ~ 30 MeV/fm 2 have 
been considered for studying the effect of quark matter nucleation on the evolution of proto- neutron stars [15| . But 
the authors in Ref. (l6| . taking into account the effects from charge screening and structured mixed phases, estimate 
70 ~ 50 — 150 MeV/fm 2 , without excluding smaller values, and an ever higher value, 70 ~ 300 MeV/fm 2 , was found 
by Afford et al. [l7j on the basis of dimensional analysis of the minimal interface between a color-flavor locked phase 
and nuclear matter. 

The surface tension for two-flavor quark matter was evaluated within the framework of the LSM by Palhares and 
Fraga [18]. In that work the authors considered the one-loop effective potential and then fitted its relevant part, which 
included both the chirally symmetric and broken state, by a quartic polynomial. The surface tension was evaluated 
using the thin-wall approximation for bubble nucleation and the estimated values cover the 5 — 15 MeV/fm 2 range, 
depending on the inclusion of vacuum and/or thermal corrections. In principle, this range makes nucleation of quark 
matter possible during the early post-bounce stage of core-collapse supernovae and it is thus a rather important result. 
It is also worth noting that a small surface tension would facilitate various structures in compact stars, including mixed 
phases fljjj . 

The present work is devoted to the evaluation of the surface tension for quark matter using both the LSM (with two 
flavors) and the NJL model (with two and three flavors) following the procedure employed in Ref. [liq . Here, the LSM 
is mainly included to check the consistency of our procedure by comparing our present results with those obtained 
by the thin- wall approximation of Ref. [18[ (we find the agreement to be very good) . The NJL model is considered 
with two and three flavors because the latter, which contains strangeness, is one of the most popular effective quark 
models used in studies related to compact stars. As explained below, the method described in Ref. [I2J makes it 
possible to express the surface tension for any subcritical temperature in terms of the free energy density for uniform 
matter in the unstable density range. Because the models employed readily provide the equation of state (EoS) for 
the full density range, they are well suited for our purpose and we may directly employ the method without any 
further approximations. In practice, the procedure is rather simple to implement and it provides an estimate for the 
surface tension that is consistent with the EoS implied by the model employed, with its specific approximations and 
parametrizations. 

The paper is organized as follows. In Sect. [TT] we review the method for extracting the surface tension from the 
equation of state. In Sec. IHII we then present the two-flavor versions of the two models considered and discuss how to 
extract the surface tension. Section ITVl is devoted to the treatment for the more realistic SU{2>) version of the NJL 
model and our numerical results are presented in Sec. [Vl both for cold matter and for temperatures up to the critical 
value. The conclusions and final remarks are presented in Sec. IVII 

II. THE GEOMETRIC APPROACH TO THE SURFACE TENSION EVALUATION 

We assume here that the material at hand, strongly interacting matter, may appear in two different phases under 
the same thermodynamic conditions of temperature T, chemical potential p, and pressure P. These two coexisting 
phases have different values of other relevant quantities, such as the energy density £, the (net baryon) density p, 
and the entropy density s. Under such circumstances, the two phases will develop a mechanically stable interface if 
placed in physical contact and it is the purpose of the present study to evaluate the associated interface tension, 7^. 

The two-phase feature appears for all temperatures below the critical value, T c . Thus, for any subcritical temper- 
ature, T < T c , hadronic matter at the (net-baryon) density pi(T) has the same chemical potential and pressure as 
quark matter at the (larger) density p2(T). As T is increased from zero to T c , the coexistence phase points (p±, T) and 
(p2, T) trace out the lower and higher branches of the phase coexistence boundary, respectively, gradually approaching 
each other and finally coinciding for T = T c . Any (p, T) phase point outside of this boundary corresponds to thermo- 
dynamically stable uniform matter, whereas uniform matter prepared with a density and temperature corresponding 
to a phase point inside the phase coexistence boundary is thermodynamically unstable and prefers to separate into two 
coexisting thermodynamically stable phases separated by a mechanically stable interface. Because such a two-phase 
configuration is in global thermodynamic equilibrium, the local values of T, p, and P remain unchanged as one moves 
from the interior of one phase through the interface region and into the interior of the partner phase, as the local 
density p increases steadily from the lower coexistence value p± to the corresponding higher coexistence value p2- 

It is convenient to work in the canonical framework in which where the control parameters are temperature and 
density. The basic thermodynamic function is thus fr{p), the free energy density as a function of the (net baryon) 
density p for the specified temperature T. The chemical potential can then be recovered as pr(p) — dpfr(p), and 
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the entropy density as st(p) — — <9t/t(p), so the energy density is £r{p) = fr{p) — Tdxfrip), while the pressure is 
P T (p) = pd p f T (p) - Mp)- 

For single-phase systems /t(p) is convex, i.e. its second derivative 9p/r(p) is positive, while the appearance of 
a concavity in /t(p) signals the occurrence of phase coexistence, at that temperature. This is easily understood 
because when /t(p) has a local concave anomaly then there exist a pair of densities, pi and P2, for which the tangents 
to fr{p) are common. Therefore fr(p) has the same slope at those two densities, so the corresponding chemical 
potentials are equal, /it(pi) = <9 p /t(pi) = c^p/xO^) = Pt{P2)- Furthermore, because a linear extrapolation of fr(p) 
leads from one of the touching points to the other, also the two pressures are equal, Pr(pi) = Pidpfripi) ~ fr(pi) — 
Pidpfo{,P?) — Jt{P2) = Pt{P2)- So uniform matter at the density p\ has the same temperature, chemical potential, and 
pressure as uniform matter at the density p 2 - The common tangent between the two coexistence points corresponds 
to the familiar Maxwell construction and shall here be denoted as f^ I {p). Obviously, /t(p) and f^ip) coincide at 
the two coexistence densities and, furthermore, /t(p) exceeds Jt{p) for intermediate densities. Therefore we have 
A/ T (p)EE/ T (p)-/^( p )>0. 

For a given (subcritical) temperature T, we now consider a configuration in which the two coexisting bulk phases 
are placed in physical contact along a planar interface. The associated equilibrium profile density is denoted by 
Pt(z) where z denotes the location in the direction normal to the interface. In the diffuse interface region, the 
corresponding local free energy density, /t(z), differs from what it would be for the corresponding Maxwell system, 
i.e. a mathematical mix of the two coexisting bulk phases with the mixing ratio adjusted to yield an average density 
equal to the local value p(z). This local deficit amounts to 

Sfr(z) = f T (z) - h - /T(P2) ' /T(Pl) (p T (z) - Pi ) , (2.1) 

Pi - Pi 

where pi is either one of the two coexistence densities. The function 5Jt{z) is smooth and it tends quickly to zero 
away from the interface where Pt{ z ) rapidly approaches pi and friz) rapidly approaches friPi)- The interface tension 
7t is the total deficit in free energy per unit area of planar interface, 

r+oo 

IT = / 5f T {z)dz . (2.2) 



As discussed in Ref. pjj], when a gradient term used to take account of finite-range effects, the tension associated 
with the interface between the two phases can be expressed without explicit knowledge about the profile functions 
but exclusively in terms of the equation of state for uniform (albeit unstable) matter, 



pp 2 (T) j 

7T = a / [2£ g Af T (p)} 1/2 ^ , (2.3) 
J P i(T) Pg 

where p g is a characteristic value of the density and £ g is a characteristic value of the energy density, while the 
parameter a is an effective interaction range related to the strength of the gradient term, C = a 2 £ g / p\ 1 . We choose the 
characteristic phase point to be in the middle of the coexistence region, p g — p c and £ g — [£o(p c ) + £ c ]/2, where £o(pc) 
is energy density at (p c ,T = 0), while £ c is energy density at the critical point (p c ,T c ). The length a as a somewhat 
adjustable parameter governing the width of the interface region and the magnitude of the tension [12j . For the LSM 
it is natural to expect that a ~ \jm a ~ 0.33 fm which, also, is approximately the value found in an application of the 
Thomas-Fermi approximation to the NJL model [20j . Therefore, we shall adopt the value a — 0.33 fm throughout the 
present work. While there is some arbitrariness in fixing these quantities, it is reassuring that the resulting surface 
tension is in excellent agreement with the value obtained in Ref. [l8| . 

With these parameters fixed, the interface tension can be calculated once the free energy density /t(p) is known 
for uniform matter in the unstable phase region, pi{T) < p < p2{T). While this is straightforward in a canonical 
formulation, where each (/?, T) characterizes only one manifestation of the system, even inside the unstable phase 
region, the task is more complicated in the commonly used grand canonical formulation because a given (p, T) phase 
point characterizes three different manifestations of the system, one stable, one metastable, and one unstable. The 
metastable solutions are located near the coexistence densities, while the unstable solutions are located in the inter- 
mediate spinodal region where uniform matter is mechanically unstable so that even infinitesimal irregularities may 
be exponentially amplified. By contrast, only irregularities of a sufficient amplitude are amplified in the metastable 
regions, leading towards either nucleation (near the lower coexistence density p{) or bubble formation (near the higher 
coexistence density p 2 )- 
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III. THE EOS FOR THE EFFECTIVE TWO FLAVOR QUARK MODELS 

In this section, we review the mean-field results for the thermodynamic potential for the two effective models when 
only two quark flavors are included. These results have been widely discussed in the literature and here we follow 
Ref. [5l[ (see Ref. [22j for results beyond the mean- field approximation). The two models are similar in the sense 
that they do not have confinement and they incorporate spontaneous chiral symmetry breaking, which happens at 
the classical level in the LSM but only via quantum corrections in the NJL model. The fermionic fields representing 
the quarks are the only degrees of freedom in the NJL model at the tree level, while the LSM also contains scalar (a) 
and pseudoscalar (tt) meson fields. 



A. The linear a model 



In standard notation, the Lagrangian density of the LSM with quarks reads 

A.SM = i(a Al 7r) 2 + i(a A1 a) 2 -C/( ( 7,7r)+^[i7^-. 9 ( ( 7 + i75r-7r)] t /', (3.1) 
where %/} is the flavor isodoublet spinor representing the quarks (u and d), and 



U{a,ir) = ^-(a 2 +TT 2 -v 2 ) 2 -Ha , (3.2) 



is the classical potential energy density. In the chiral limit (obtained for H — 0) the chiral symmetry, SU(2)v x SU(2)a, 
is spontaneously broken at the classical level, and the pion is the associated massless Goldstone boson. For H ^ 0, 
the chiral symmetry is explicitly broken by the last term in U (cr, 7r) which gives the pion a finite mass at vanishing 
T and p. The parameters are usually chosen so that chiral symmetry is spontaneously broken in the vacuum and 
the expectation values of the meson fields are (cr) = /„■ and (tt) = where f n = 93 MeV is the pion decay constant. 
Following Ref. [U, we fix the parameters as follows: v 2 ~ (87.73 MeV) 2 , A 2 ~ 20, and H ~ (12.1 GeV) 3 . Using the 
standard relations, H = f^m 2 , v 2 = f 2 — m 2 /A 2 , and m 2 = 2A 2 / 2 , we obtain the meson masses, = 138 MeV and 
m a = 600 MeV. The coupling constant g is usually fixed so that the effective quark mass in vacuum, M vac = gf„, 
be about one third of the nucleon mass, which gives g ~ 3.3. We note that the same parameter set was also used 
to evaluate the surface tension, jx> i n Ref- flal- T o the one- loop level, the grand canonical potential is obtained by 
integrating the action over the fermionic fields [21j . 

n LSM (a,TT;T :f i q ) = U(a,TT)-2N f N c J-0^{E-Tln[l-?i+]-T\n[l-n-]} , (3.3) 

where iV c = 3, Nt — 2, E 2 —p 2 + AI 2 , and n = {1 + exp[(E =p p^/T]}^ 1 represent the particle/antiparticle thermal 
occupancies with p q = p u = pa, where p = 3p q . For for given values of T and p, the equilibrium values of the meson 
fields are obtained by minimizing f2(cr, tt; T, p q ) with respect to those, yielding the most likely values a and jr. The 
latter one vanishes in the mean-field approximation, so the associated constituent quark mass is given as M 2 — g 2 a 2 . 
The minimum value of the grand potential represents minus the equilibrium pressure, f2 m i n (T, fi) = —P, so the net 
quark density is given by p q = (dP/dfi q )T, where p q = 3p. The entropy density given by s = (dP/dT)^ , while the 
energy density, £ , can then be obtained by means of the standard thermodynamic relation P = Ts — £ + fip and the 
free energy density is / = £ — Ts = up — P. 

In the neighborhood of the phase coexistence line in the (p, T) plane, the grand potential has three extrema 
representing stable, metastable, and spinodally unstable matter. As emphasized above, the extraction of the surface 
tension by the geometric approach requires the consideration of all three extrema. 

In contrast to the NJL model, the vacuum term represented by the first term in the integrand of Eq. (I3.3[) is not 
essential for the spontaneous breaking of chiral symmetry. In the LSM this already happens at the classical level and 
the symmetry restoration is driven mainly by the terms containing n^ 1 . Therefore, we neglect the vacuum term in 
the present LSM application, where the aim is to compare our estimates with the zero-temperature interface tension 
obtained in Ref. [18|], 70 — 12.98 MeV/fm 2 , where the relevant part of the same thermodynamical potential was fitted 
with a quartic polynomial. In our approach such a fitting procedure is not necessary because the thermodynamic 
potential is evaluated for all values of p and T. This will lead to somewhat different numerical values for the surface 
tension 77-. It was shown in Ref. [l8| that the inclusion of vacuum terms at T = increases the surface tension value 
from 70 = 12.98 MeV/fm 2 to about 70 ~ 17MeV/fm 2 . In practice, further refinements including vacuum and in- 
medium two-loop corrections are possible by following the same technical steps that were employed in the evaluation 
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It is now straightforward to determine the phase-coexistence line in the (//, T) plane which forms the starting point 
for determining all quantities related to Jt- It starts at f/i = 918 MeV, T — 0) and terminates at the critical point 
(fj, c = 621 MeV, T c — 99 MeV), which agrees with Ref. '21}. We now have all the ingredients needed for determining 
the coexistence densities pi(T) and P2(T) as well as the characteristic values p g and £ g appearing in Eq. (|2.3[) . 

The difference A/(p) can be readily determined numerically by considering the stable (global) minimum, the 
metastable (local) minimum, and the unstable (local) maximum appearing in the thermodynamical potential, as will 
be explicitely shown in section |Vj Although our method for extracting the surface tension does not require the profiles 
a(z) and p{z), these functions do provide interesting additional information about the interface. To obtain the profile 
functions within the LSM, it suffices to consider the grand canonical potential in the a direction only, i.e. taking 
tv = 0. At T = it can be expressed in terms of the Fermi momentum pp (given by p 2 F = p? — g 2 o ), 



«LSM(a, 7T = 0; T = 0, p) = U(a) - \ 2p qV % Ham r 



PqPF - (gcr) 2 In 



Pf + Pg 

go- 



(3.4) 



We now employ the local density approximation, so the local Fermi momentum, pp(r), is related to the local density, 
p q (r) = (ip + ip), by p q = (N c Nf /3ir 2 )p F (r). Furthermore, the local scalar density p s (r) = (i/iip) is given by 



(3.5) 



where we have used that the local Fermi energy Ep(r) — y/pp(r) 2 + M(r) 2 equals the (constant) chemical potential 
p q = p/3. We note that 

[dMLSu(a,O;0,p)l =a{r) - U'(a(r)) + gp s (r) . (3.6) 

Then the stationary Euler-Lagrange equation with tt = and iptp replaced by p s (^), provides an equation for the 
local value of the order parameter, er(r), 

V 2 a(r)~[U'(a(r))+gp s (r)} = 0. (3.7) 

For semi-infinite geometry, the profile of the order parameter, cr(z), can then be obtained by solving the corresponding 
Euler-Lagrange equation, 

dla(z) - [c^lsmMU^) = , (3.8) 

with the boundary conditions that the order parameter approach the zero-temperature coexistence values far from 
the surface, 

a(z -> -co) -> o- 1 (T = 0) = / 7r , cr(z -» +oo) -> a 2 (T = 0) = 0.13 f„. (3.9) 

Once (j{z) is known, so is the mass M(z), and we can then obtain the local Fermi momentum pf{z) and, conse- 
quently, the local net baryon density p(z). One may then define the associated interface location function |24j|. 

/ \ _ d z p(z) 

g{z) = , (3.10) 

P2 - Pi 

which allows us to obtain additional information, such as the mean interface location, z = (z) = J zg(z)dz, the 
associated interface width b, where b 2 = ((z — z) 2 ), as well as a measure of the profile skewness which is given by the 
dimensionless parameter 73 = ((z — z) 3 )/b 3 [HJ|. The calculated results for cr(z), p(z), and g(z) are shown in Fig. [T] 
The origin is conveniently located at (z), the width is b = 1.85a (where a = l/m a m 0.33 fm) and the skewness is 
73 = 0.6 . As the temperature is increased, the profiles widen progressively and grow more symmetric, as also found 
in Ref. El. 



B. The Nambu— Jona-Lasino model 

We now consider the standard version of the two-flavor NJL model [25|. Its Lagrangian density is based on a 
chirally symmetric four-fermion interaction, 



£njl = i/j {il^ —m)ip + G [{Uf + {i>i l5 Ti>) 2 ] , 



(3.11) 
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FIG. 1: The order parameter a(z) in units of (dashed), the net baryon density p(z) in units of po (solid), and the surface 
location function g(z) (dotted), as functions of the dimensionless depth (z — (z))m a relative to the average surface location. 



where ip is to be interpreted as in the LSM. Furthermore, it is assumed that m u = wid so the mass matrix is given by 
m c = mdiag(l, 1). In the mean-field approximation, the grand canonical potential reads (2ll l22l [26j 

^n.tl(m,T) = (M ~ G m)2 ~ WfNc J {E - Tin [l - n+] - Tin [l - n~] } , (3.12) 

with the same definitions as used in the LSM. For each value of T and p the dynamical mass is of the form Mdiag(l, 1) 
because of the assumption of isospin symmetry (m u — md — m ) an d chemical equilibrium (p u = pd = p); alternative 
scenarios may also be considered [27] . The single dynamical mass M is then obtained by minimizing f2 with respect 
to M, leading to the well known gap equation, 

M = m + 2GN c N f j^^[l-n--n-] . (3.13) 

Although the thermodynamic potentials for LSM (Eq. ()3.3|) ) and the NJL model (Eq. l3.12"T) ) have the same structure 
as far as the loop contribution is concerned, some important differences between the two models exist. First, we note 
that within the NJL the quark mass acquires its constituent value only when quantum corrections (loop terms) are 
computed. Therefore, contrary to the LSM, the divergent term represented by the second term on the right-hand side 
of Eq. (|3.12[) plays a central role regarding the (dynamical) chiral symmetry breaking. Another difference between 
the LSM and the NJL model, in 3+1 dimensions, is that the latter is not renormalizable since the coupling G carries 
dimensions (energy" 2 ). This means that potential divergencies cannot be systematically eliminated by a redefinition 
of the original parameters. By considering it as an effective model, one gives up the very high energies and evaluates 
all the integrals up to an ultraviolet (non-covariant) cutoff A, as the notation in Eqs. (|3.12[) and (|3.13[) implies. Then 
A is treated as a "parameter" which will be fixed, together with G and to, so as to yield the values of physical 
observables such as m^, f n that reproduce the phenomenological value of (ifiip). For example, in Ref. [2l[ the authors 
reproduce f„ = 93 MeV and = 138 MeV using A = 631 MeV and GA 2 = 2.19 with to = 5.5 MeV. These parameter 
values, which we label "set I", predict a first-order phase transition starting at T = 0,/i = 1045.5 MeV and ending 
at the critical point (T c = 46 MeV, p c = 996 MeV), while the constituent quark mass in vacuum is M vac = 337 MeV. 
In their study of the chiral phase transition in the presence of spinodal decomposition the authors in Ref. (2(| use 
A = 587.9 MeV and GA 2 = 2.44 with m = 5.6 MeV in order to reproduce f v = 92 A MeV and = 135 MeV obtaining 
M vac = 400 MeV. We shall also consider these parameter values, which we label "set II", in order to estimate the 
influence of different parametrizations in the estimation of jt- With parameter set II the first-order transition line 
starts at (T = 0,/x = 1146.3 MeV) and ends at (T c = 81 MeV, ^ c = 990 MeV). In general, a larger value of GA 2 
enlarges the coexistence region. As in the LSM case, the quantities pi, p 2 , p g , £ g , and A/(p) entering the expression 
(|2.3[) for the surface tension can be obtained from the equation of state. As already emphasized, the numerical value 
for the length scale a is chosen to be l/m a ~ 0.33 fm (which is about the value found in a Thomas-Fermi application 
to the NJL model [20]). The remaining two numerical inputs, p g and e g , are automatically fixed once the EoS has 
been determined. 
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FIG. 2: The constituent quark masses as functions of of the baryon chemical potential \x for T = in the three-flavor NJL 
model. The top curve shows M s , while the curve represents M u = Md- Phase coexistence occurs at /i = 1083 MeV and the 
corresponding mass values are indicated by the dots. These are are joined by the Maxwell lines (dashed) that trace out the 
gradual phase conversion when full equilibrium is maintained; M 3 then decreases from 549 to 464.4 MeV, while M u = Md 
decreases from 367.6 MeV to 52.5 MeV. Between the two dotted lines there are three solutions for a given \x: stable matter, 
metastable matter, and (between the two extrema) spinodally unstable matter. 
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FIG. 3: The pressure, P, as a function of the compression p/po for the NJL model as obtained with the two parameter sets. 
The two coexistence points (dots) are joined by the Maxwell line (dashed) along which the global equilibrium evolves as the 
density is increased through the phase coexistence region. 



IV. THE EOS FOR THE NJL MODEL WITH THREE QUARK FLAVORS 

In stellar modeling, the structure of the star depends on the assumed EoS built with jappropriate models while the 
true ground state of matter remains a source of speculation. It has been argued [28h31| that strange quark matter 
(SQM) is the true ground state of all matter and this hypothesis is known as the Bodmer-Witten conjecture. Hence, 
the interior of neutron stars should be composed predominantly of u, d, s quarks (plus leptons if one wants to ensure 
charge neutrality which is not the case in the present work). Strangeness is implemented in the SU(3) version of the 
NJL model which is given by 

8 

C = j(i lfl d» - m)iP + GJ2 [WA q ^) 2 + (^7sA» 2 ] -K {det / [Vi(l + 7 5 )V]+det / [Vi(l-7 5 )V]} , (4.1) 

a=0 

where "0 = ( u ,d,s) T denotes a quark field with three flavors (and three colors), and m = diag^ (m u , mj, m s ) is the 
corresponding mass matrix. Here we assume m u = md ^ m s indicating that isospin symmetry is observed while the 
SU(3) flavor symmetry is explicitely broken. The eight Gell-Mann matrices are represented by A a (a = 1, ...,8) and 
A = y/2/3I. More details concering this version of the NJL model can be found in Ref. [27|. In the mean-field 
approximation the thermodynamical potential is given by 



n N j L3 (T,// g ) = [^m 5 (T,m 9 ) + 2G^] - AK<t> u <t> d <t> s 

q—u,d,s 



(4.2) 
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FIG. 4: The baryon chemical potential, p, as a function of the compression p/po for the NJL model as obtained with the two 
parameter sets. The two coexistence points (dots) are joined by the Maxwell line (dashed) along which the global equilibrium 
evolves through the phase coexistence region. 
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FIG. 5: Left panel: The free energy density at zero temperature fo(p) (solid) and the associated Maxwell line fo I {p) (dashed) 
as functions of the compression p/po for the NJL model with parameter set II. The Maxwell line is tangent to /o(p) at the two 
coexistence densities pi = and p2 = 3.027 po- Right panel: The quantity A/o(p) as a function of p/po obtained with the NJL 
model for parameter set I (continuous) and parameter set II (dot-dashed). 



The term ^Im 9 , which represents the contribution of a gas of quasiparticles with mass M q , is given by 

n Mq = -2N C j £L { Eq - Tin [1 - n+] - Tin [l - n~] } , (4.3) 
Jp<\ { Z7T J 

where E 2 = p 2 + M 2 and rv^ = {1 + exp[(E q =p /i g )/T]} _1 represent the particle/antiparticle distribution function. 
For the quark condensates, 4> q = {ipqipq}, one has 

Finally, the gap equation is 



Mi = rrii — 4G(j)i +2K(f>j(pk , {i, j, k) = any permutation of (u,d,s) , (4-5) 

which contains a non- flavor mixing term proportional to G as well as a flavor mixing term proportional to K . In our 
numerical analysis we adopt the parameter values of Ref. [32[ which are m u — md — 5.5 MeV, m s — 140.7 MeV, GA 2 = 
1.835, KA 5 = 12.36, and A = 602.3 MeV. Then, at T = and fi f = 0, one reproduces f n = 92.4 MeV, = 135 MeV, 
rriK = 497.7 MeV, and m n i = 960.8 MeV. For the quark condensates one obtains (j) u = <t>d — —(241.9 MeV) 3 , and 
4> s = -(257.7MeV) 3 . The constituent quark masses are then given by M u = M d = 367.7MeV and M s = 549.5 MeV. 
The pressure, P, and energy density, £, follow from the usual expressions, 

P = -0 N jl 3 (T, K» and P = Ts-£+ J2 NPi ■ ( 4 -6) 

q— u,d,s 
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FIG. 6: The surface tension 7t (relative to its zero-temperature value 70) as a function of the temperature T (measured relative 
its critical value T c ), as obtained with the NJL model for parameter set II. 

Here, for simplicity, we take the chemical equilibrium condition, p u — pd — p s = p q = p/3 which yields M u = M<i 
also at finite T and/or p. Of course, for a realistic description of neutron-star matter, charge and strangeness neutrality 
need to be taken into account, which is technically straightforward. 

The phase diagram for this three-flavor quark model model in the p-T and p-T planes can be found in Ref. 
[33| : the first-order transition line starts at (T = 0, p, = 1083 MeV) and ends at the critical point (T c = 67.7 MeV, 
p c = 955.2 MeV). As already discussed, we need the EoS inside the phase coexistence region which can be obtained 
by examining how the effective masses behave in this domain. This behavior is shown in Fig. [2] for T — 0; these 
results go beyond those of Ref. [26| by also considering the strange quark mass. To understand this figure, let us 
recall that, in most situations, one is generally interested only in those solutions of the gap equation that correspond 
to global (stable) minima of the thermodynamical potential. However, when a first-order phase transition is present, 
there are two different such solutions for the same thermodynamic conditions of temperature, chemical potential, and 
pressure (corresponding to the solid dots on Fig. [5]). As the net baryon density (which serves as a convenient order 
parameter) is increased from its lower coexistence value p\ to its higher coexistence value P2, the thermodynamically 
favored state is a Maxwell mixture of the two coexisting phases and the overall average of the energy per net baryon 
or the effective mass, for example, evolve monotonically along the socalled Maxwell line, as the composition of the 
mixture changes from being entirely one phase to being entirely the other. In the region between the dotted lines 
the gap equation has three solutions, leading to the back-bending evolution brought in that diagram. It is precisely 
this typical first-order behavior that will be reflected in the thermodynamical quantities, such as the pressure and 
densities, as Figs. [3] and [4] show. This behavior is responsible for the fact that there is a (positive) deviation A/(p), 
which then in turn leads to the surface tension. 



V. NUMERICAL RESULTS 

We now turn to our numerical results for the surface tension jt- To this end we need to determine the free energy 
density fr(p), which requires the evaluation of Pt{p) and pt{p) for uniform matter thermodynamically unstable 
region of the phase diagram. For the considered temperature T, the associated density region is bounded by the two 
coexistence densities p\ and P2, for which the chemical potential p has the same value, as does the pressure P. As the 
density p is increased through the lower mechanically metastable (nucleation) region, p and P rise steadily until the 
lower spinodal boundary has been reached. Then, as p moves through the mechanically unstable (spinodal) region, 
both p and P decrease until the higher spinodal boundary is reached. They then increase again as p moves through 
the higher mechanically metastable (bubble- formation) region, until they finally regain their original values at p = p2- 
It is convenient to express the (net) baryonic density p in units of the nuclear saturation density, po = 0.153/fm 3 . 
Generally, as is common practice, we subtract from the pressure any finite value it may have in the vaccum. 

A. Zero temperature 

Let us start with T = for which the relevant results can be readily obtained by taking the T — > limit in Eqs. 
(|3.3p and (13.12j) (see, e.g. Refs. [H, H3, H3 ) • Figure [3] shows the pressure as a function of the degree of compression 
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p/po obtained with the NJL model for both parameter sets I and II; the latter has a stronger coupling and a larger 
coexistence region. A qualitatively similar behavior is observed in Fig. 3] which shows p as a function of p/ po for the 
NJL model with both parameter sets. Figure [5] shows the behavior of the free energy fo(p) and its corresponding 
Maxwell line /q (p) for parameter set II. In the right panel of Fig. [5] we display the difference between these two free 
energies, A/(p) = fo(p) — fn 4 {p) for both parameter sets; this is the key quantity for the determination of the surface 
tension. The LSM and the three-flavor NJL model yield similar results. At temperatures below criticality, T < T c , 
the thermodynamical potential has two degenerate minima determining the densities of the two coexisting phases, 
pi and pi. In all cases studied here, the lower coexistence density vanishes, p\ = 0. As for the higher coexistence 
density, the LSM yields P2/P0 — 1-54, the two-flavor NJL model yields pi/po — 2.13 with set I and P2/ Pa — 3.03 
with set II, while the three-flavor NJL model gives P2/P0 — 2.62. It clear from this figure that set I produces a much 
weaker phase transition because A/(p) is much smaller that for set II, as is indeed reflected in the jt values shown 
in Table I. In fact, set II produces a greater coexistence region (pi ~ 0, P2 — 3.03,po,T c — 81MeV) when compared 
to set I (pi ~ 0, P2 — 2.13po, T c — 46 MeV), which is in accordance with the well known fact that set II should cause 
the size of the first-order transition line to be longer than the one produced by set I which has a weaker coupling. 
Further refinements, such as finite- N c corrections [371 ] . contributions from thermal flucutations, and the inclusion of a 
repulsive vector interaction [38] , also tend to shrink the first-order transition line [22] so that, within a fixed parameter 
set, one should expect these effects to reduce jt- 

Table I summarizes all our results for 70 and also lists the characteristic values £ g and p g as well as the location 
of the critical point (T c , p c ). The table also provides information related to thermodynamic potential at T = and 
p = by showing the values of the constituent quark mass in vacuum (M vac ), which is related to the distance from 
the global minimum to the origin, as well as the bag constant which gives the energy difference between the local 
maximum and the global minimum of the potential in vacuum; these values were taken from Refs. [2TL [27j . 



Model 


70 


Bo M™ d c 


M vac 


T c p c 


Pg/po 


£* 


NJL (I) 


7.11 


100 337 




46 996 


2.00 


342.85 


NJL (II) 


30.25 


141.4 400 




81 990 


2.42 


495 


LSM 


13.18 


60 306.9 




99 621 


1.19 


219.25 


NJL3 


20.42 


291.7 367.6 


549.5 


67.7 955.2 


1.87 


326.8 



TABLE I: Summary of inputs and results. The length parameter was taken as a = 0.33 fm. The zero-temperature bag constant 
Bo and the characteristic energy density £ g are given in MeV/fm 3 , while the effective quark masses in vacuum Mj Vac as well 
as the critical values (i c and T c are in MeV. The resulting zero-temperature surface tension 70 (given in MeV/fm 2 ) may be 
compared with the value 70 = 12.96 MeV/fm 2 obtained in Ref. [13]) with the LSM in the thin-wall approximation. 

We finally note that Palhares et al. (lSf . using the approximation jt ~ / \d z cr(z)\ 2 dz obtained the estimate 
70 ~ 12.98 MeV/fm 2 which is very close to our LSM value of 13.18 MeV/fm 2 and also rather similar to the value 
12.19 MeV/fm 2 resulting from evaluating that intergal using our LSM profile function cr (z). 

B. Finite temperature 

One can easily consider finite temperatures within the employed models. The interface tension is expected to 
decrease with increasing temperature because both the coexistence densities pi and the associated free energy densities 
fr(pi) move closer together at higher T; they ultimately coincide at T c where, therefore, the tension vanishes. This 
general behavior is confirmed by our calculations, as shown in Fig. [6] The LSM, the NJL with parameter set I, and the 
three-flavor NJL display similar behaviors. The temperature dependence of the surface tension may be relevant for the 
thermal formation of quark droplets in cold hadronic matter found in "hot" protoneutron stars whose temperatures, 
T*, are of the order 10-20 MeV [H, H3, l35j ]. For T* the relevant value of jt, may be estimated by using table I 
together with Fig. [SJ For example, the three-flavor NJL model yields 7t» ~ 14 — 18 MeV/fm 2 . The temperature 
dependence of the surface tension is also important in the context of heavy-ion collisions, because it determines the 
favored size of the clumping caused by the action of spinodal instabilities as the expanding matter traverses the 
unstable phase-coexistence region. 

VI. CONCLUSIONS 

In this work we have shown that the interface tension related to a first-order phase transition may be evaluated 
once the uniform-matter equation of state is available for the unstable regions of the phase diagram. It is a convenient 
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feature of the method employed that knowledge of the interface profile functions is is not required, because their 
determination can be quite complicated, as is the case for NJL model [2(| (although it is easy for the LSM). In 
addition to the EoS, the geometrical approach also requires a proper setting of three input parameters, namely the 
characteristic densities p g and £ g together with the length scale a. While this does encumber the numerical results 
with some degree of uncertainty, our zero-temperature LSM result, 70 = 13.18 MeV/fm 2 , agrees within a few percent 
with the approximate value obtained in Ref. [18|, thus suggesting that those parameters were chosen reasonably. 

The surface tension determined in the present fashion is entirely consistent with the employed model, including 
the approximations and parametrizations adopted. For the effective quark models employed here, this amounts to 
considering all the solutions to the gap equation (stable, metastable and unstable) and determine the relevant effective 
quark masses. In most non-perturbative approximations (large N c , mean field, etc.) the various quantities of interest, 
such as the free energy density, become functions of this effective mass and will therefore also reflect the metastable 
and unstable character of the configuration considered. As a cross check on our procedure, we have evaluated 70 for 
the LSM obtaining a result that differs by only about 2% from estimates based on the thin- wall approximation fl8| . 
We have investigated the two-flavor NJL model as well as its more realistic three-flavor version. 

Our main conclusion is that all these effective models generate relatively low values for the the surface tension. This 
would favor the formation of quark matter and may thus have important astrophysical consequences regarding the 
existence of pure quark stars. Of particular interest is the three-flavor NJL result, 70 = 20.34 MeV/fm 2 , because this 
model is widely used in studies related to neutron stars. Here, for simplicity, we have considered pure quark matter 
where all flavors share the same chemical potential, but it is just a technical matter to generalize our procedure so 
as to include leptons (e, fi) in order to enforce /3 equilibrium (fi e i = fi s = jj, u + // e , /j, e = /z^), although the additional 
chemical potential introduces an increased degree of complexity into the features of the phase transition. 

In principle, more refined treatments, such as the Polyakov-NJL model, can also be considered within the same 
framework. However, because the effects of the Polyakov loop become more important above 100 MeV [36[ we 
believe that our results, especially the three-flavor ones, can be considered as reasonably accurate, although numerical 
variations may arise due to the parametrizations and approximations adopted. 
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